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ABSTRACT 


Ionospheric  distortions  limit  the  usable  signal  bandwidth  of  HF 
over-the-horizon  paths.  By  measuring  the  transfer  function  of  the 
path  and  correcting  for  it  in  real  time  improved  bandwidth  capability 
results.  To  determine  the  feasibility  of  such  a  real-time  correction 
technique,  data  has  been  gathered  on  an  HF  link  operated  by  Stanford 
University  and  analyzed  by  computer  at  MITRE.  A  non-real  time  computer 
simulation  of  a  correction  technique  has  shown  that  for  the  bandwidths 
analyzed,  the  correction  technique  is  feasible  and  that  the  corrections 
will  not  deteriorate  significantly  for  several  seconds. 
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SECTION  I 


INTRODUCTION 

Transmission  of  HF  signals  via  ionospheric  refraction  is  limited 
by  ionospheric  distortions.  These  distortions  appear  as  dispersion, 
absorption  fading,  multiple  ray  paths  and  faraday  rotation,  all  of 
which  vary  with  time.  Compensation  for  distortion  would  permit  the 
use  of  increased  bandwidth,  thereby  effecting  multipath  isolation 
and  improving  resolution,  accuracy  and  clutter  rejection.  The  aim  of 
MITRE  adaptive  processing  techniques  program  has  been  to  measure  the 
characteristics  of  this  medium,  devise  techniques  that  are  capable 
of  correcting  for  the  distortions,  and  build  equipment  which  tracks 
or  adapts  to  changes  in  the  ionospheric  characteristic  in  real  time. 

Early  in  the  program  a  theoretical  adaptive  processing  technique 
was  devised.  Simultaneously  it  was  established  that  suitable  ionos¬ 
pheric  data,  to  which  the  adaptive  processing  technique  could  be  applied, 
was  not  available  but  that  an  ionospheric  sounder  capable  of  producing 
the  required  data  was  operational  between  Lubbock,  Texas  and  the 
Stanford  Electronics  Laboratory  (SEL)  of  Stanford  University.  This 
sounder  was  ideally  suited  because  of  its  slowly  sweeping,  phase 
coherent,  linear  FM  waveform.  Since  MITRE  possessed  digital  recording 
equipment  and  a  library  of  digital  signal  processing  programs,  a 
cooperative  experimental  venture  between  MITRE  and  Stanford  was  entered 
into.  Data  gathering  exercises  were  performed  in  February  and  June 
of  1967. 
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This  paper  presents  distortion  measurements  made  over  this  one- 
hop  oblique  path  and  the  results  of  a  computer  simulation  of  a 
correction  technique  as  applied  to  these  data. 
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SECTION  II 


THEORETICAL  ADAPTIVE  PROCESSING  TECHNIQUES 

If  the  ionosphere  is  modeled  as  a  linear  time  invariant  filter 
(time  invariant  over  the  time  required  for  measurement  and  correction), 
its  characteristic  can  be  represented  by  a  transfer  function 


1(f)  =  E(f )  e 


jD(f) 


(1) 


where  E(f)  is  an  arbitrary  amplitude  characteristic  as  a  function  of 
frequency  and  D(f)  is  an  arbitrary  phase  characteristic  as  a  function 
of  frequency. 

The  inverse  to  the  ionospheric  transfer  function,  if  it  is 
realizable,  can  be  represented  as  the  reciprocal  of  the  amplitude 
characteristic  and  the  complex  conjugate  of  the  phase  characteristic, 
as  shown  in  Equation  2. 


I. 

inverse 


(o 


E(f )  6 


(2) 


This  inverse  function  is  the  basis  for  the  corrections  described  in 
this  paper.  A  receiver  processor  capable  of  assuming  this  inverse 
transfer  function  would  be  able  to  correct  for  all  undesired  signal 
distortions.  The  mismatch,  subsequent  degradation  in  signal- to-noise 
ratio  and  parameter  distortion  associated  with  any  inverse  filter  in 
the  presence  of  multipath  would  not  be  acceptable  for  all  applications. 
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In  these  instances  compensation  for  dispersion,  but  not  multipath, 


can  be  accomplished  by  applying  a  correction  derived  from 
fit  to  the  phase  characteristic.  (This  technique  has  the 
of  approximating  matched  filter  conditions.) 


a  smoothed 
advantage 
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SECTION  III 


IONOSPHERIC  TRANSFER  FUNCTION  MEASUREMENT  TECHNIQUE 

The  ionospheric  transfer  function  could  be  obtained  from  the 
ionospheric  impulse  response.  However,  it  is  difficult  to  contain 
enough  energy  in  a  narrow  impulse.  Or  the  ionospheric  transfer  function 
could  be  obtained  from  the  amplitude  and  phase  of  CW  frequencies  spaced 
across  the  frequency  band  of  interest.  The  CW  frequencies,  however, 
must  all  be  phase  and  amplitude  calibrated. 

The  technique  of  this  experiment  is  made  possible  because  of  the 
development  by  Hewlett  Packard  and  Stanford  Electronics  Laboratory  of 
a  phase  coherent  frequency  synthesizer. ^  This  synthesizer  can  be 
programmed  to  step  at  a  very  fast  rate  from  one  frequency  to  the  next 
in  small  frequency  steps  and  in  a  phase  coherent  manner.  An  excellent 
approximation  to  a  linear  frequency  modulation  is  thereby  generated. 
This  linear  FM  signal  is  transmitted  through  the  ionosphere  and  upon 
reception  is  demodulated  by  mixing  it  with  a  delayed  replica  of  the 
transmitted  waveform.  Normally,  such  a  correlation  process  yields  a 
constant  frequency  directly  related  to  the  time  delay  from  transmitter 
to  receiver.  However,  due  to  ionospheric  distortion  the  output  signal 
contains  amplitude  and  phase  modulation.  In  fact,  due  to  the  linear 
time -frequency  relation  of  the  sounding  signal,  the  amplitude  and  phase 

^HF  Measurement  Using  Extended  Chirp-Radar  Techniques,  by  R.  B.  Fenwick 
and  G.  H.  Barry.  Stanford  Electronics  Laboratory,  Technical  Report 
No.  103,  June  1965. 
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modulation  which  emerges  from  the  correlation  mixer,  as  a  function  of 
time,  is  the  amplitude  and  phase  transfer  function  (Equation  1)  of  the 
ionosphere.  Thus,  a  direct  measure  of  the  ionospheric  transfer  function 
can  be  obtained  by  taking  quadrature  samples  of  the  output  waveform  at 
equal  time  intervals.  This  is  true,  because  for  large  time-bandwidth 
product  linear  FM  waveforms  the  amplitude  is  nearly  constant  and  the 
phase  is  nearly  quadratic  as  a  function  of  either  time  or  frequency. 

When  the  ionospheric  distortions  are  included  in  the  path,  it  is  only 
necessary  that  the  phase  and  amplitude  distortions  be  slowly  varying 
with  respect  to  the  phase  variation  of  the  modulating  waveform. 

There  are  three  major  advantages  of  the  linear  FM  waveform.  Phase 
calibration  is  automatic  once  the  linear  FM  waveform  has  been  achieved. 
Energy  is  transmitted  continuously  providing  high  energy  per  sample, 
and  narrow  band  signals  present  in  the  HF  band  are  discriminated 
against. 
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SECTION  IV 


ADAPTIVE  PROCESSOR  REALIZATION 

The  manner  in  which  the  above  described  measurement  can  be  applied 
to  correct  or  compensate  for  distortion  depends  heavily  upon  the  modu¬ 
lation  that  is  being  transmitted  and  received.  For  example,  to  have 
a  real-time  processor  capable  of  receiving  arbitrary  waveforms 
requires  the  implementation  and  control  of  an  adaptive  filter  whose 
transfer  function  is  the  inverse  of  the  ionospheric  path.  In  this 
general  case,  a  tapped  delay  line  with  amplitude  and  phase  weights 
on  each  tap  could  be  continuously  controlled  to  provide  the  inverse 
transfer  function.  Or,  a  bank  of  band-pass  filters  with  amplitude  and 
phase  weights  on  each  filter  could  be  employed. 

Although  linear  FM  is  an  ideal  waveform  to  sound  the  ionosphere, 
it  would  not  be  necessary  to  use  this  waveform  for  subsequent  trans¬ 
missions  if  the  generalized  inverse  filter  were  implemented.  However, 
if  the  signal  transmitted  consists  of  a  frequency  vs  time  function  and 
if  the  receiver  consists  of  a  correlation  mixer,  a  programmed  oscillator, 
and  a  spectrum  analyzer,  then  a  more  easily  implemented  alternative  to 
the  generalized  inverse  filter  can  be  employed.  In  such  an  instance, 
the  phase  correction  can  be  made  by  modifying  the  receiver  local 
oscillator,  and  the  amplitude  correction  can  be  made  by  a  time-varying 
gain  control  in  the  receiver  channel.  Since  the  initial  goal  of  the 
MITRE  program  is  to  demonstrate  the  feasibility  of  ionospheric  distortion 
correction,  the  technique  of  active  control  of  receiver  gain  and  phase 
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is  employed.  Linear  FM  is  used  both  for  measuring  the  ionospheric 
path  and  for  subsequent  transmissions  to  be  corrected.  The  received 
data  is  recorded  and  the  correction  is  applied  in  a  computer  in  non- 
real  time. 
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SECTION  V 


MITRE -SEL  EXPERIMENT  AND  DATA  ANALYSIS 

Figure  1  shows  a  block  diagram  representation  of  the  joint 
MITRE-SEL  experiment.  The  Texas  transmitter,  the  ionosphere,  and  the 
Stanford  demodulator  are  shown.  The  processing  and  recording  equip¬ 
ment  built  by  MITRE  is  shown  within  the  dotted  lines.  The  IF  signal 
is  mixed  in  quadrature  mixers  to  approximately  DC.  After  low  pass 
filtering,  these  signals  are  sampled,  converted  to  digital  form  and 
stored  on  magnetic  tape.  Separating  the  signal  into  quadrature  sine 
and  cosine  components  allows  both  amplitude  and  phase  to  be  sampled 
at  a  single  point  in  time.  Once  the  data  is  in  digital  form,  further 
equipment  distortions  are  avoided  and  manipulations  of  any  type  can 
be  performed  in  a  digital  computer. 

The  particular  implementation  employed  a  digital  stepping  recorder 
whose  normal  data  storage  rate  limited  the  bandwidth  of  the  received 
IF  signal  to  50  Hz.  The  signal  frequency  band  of  interest  was 
positioned  within  the  recorder  bandwidth  by  selecting  an  appropriate 
local  oscillator  frequency.  The  signal  recorded  was  usually  received 
one-hop  via  the  F  layer. 

Data  were  recorded  over  a  period  of  several  hours  on  February  25, 
1967.  During  this  session  nine  frequency  regions,  each  350  kHz  wide 
were  swept  repetitively  for  twelve  seconds.  Of  the  34  records  analyzed 
two  examples  are  included  here.  The  sweep  rate  was  1  MHz/ second. 

Figure  2  shows  the  received  signal  for  the  first  of  a  set  of  sweeps 
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Figure  I.  SEL/MITRE  OBLIQUE  SOUNDING  EXPERIMENT 
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Figure  2.  AMPLITUDE,  PHASE  AND  PHASE  RESIDUALS  vs  FREQUENCY 


between  20.15  MHz  and  20.5  MHz.  The  maximum  usable  frequency  (MUF) 

was  approximately  26  MHz  at  this  time.  Figure  2a  shows  amplitude  and 

Figure  2b  shows  phase  versus  frequency  for  the  same  sweep,  while 

Figure  2c  shows  the  phase  residuals  after  a  best  straight  line  (second- 

order  polynominal)  fit.  The  amplitude  and  phase  residuals  shown  are 

a  measurement  of  the  ionospheric  transfer  function.  The  amplitude 

recording  shows  a  fluctuation  of  about  six  cycles  in  the  band  swept 

as  well  as  a  lower  perturbation  rate  of  about  1.3  cycles.  The  maximum 

and  minimum  amplitude  values  differ  by  about  4  to  1.  The  phase 

residuals  curve  roughly  approximates  a  parabola  with  a  maximum  phase 

excursion  of  259  degrees.  A  perfectly  parabolic  phase  residual  would 

imply  a  perfectly  linear  delay  dispersion.  For  this  assumption,  the 

delay  dispersion  can  be  computed  from  the  phase  distortion  knowing 

the  linear  FM  slope  of  1  MHz  per  second.  The  dispersion  and  the 

maximum  phase  excursion  0  are  related  as  follows: 

max 

dispersion  -  360  (^dvldth)  <3) 

where  bandwidth  is  350  kHz  in  this  case,  and  0  is  in  degrees. 

Solving  Equation  (3)  for  this  example,  the  dispersion  is  equiva¬ 
lent  to  16.5  ^isec  of  delay.  The  Fourier  transform  corresponds  to  a 
time  domain  impulse  response.  Figure  3  shows  the  transform  of  this 
sweep.  A  40  db  Taylor  weighting  was  used  to  suppress  sidelobes. 

The  16.5  [-lsec  of  dispersion  calculated  assuming  a  perfectly  parabolic 


12 


Figure  3.  WEIGHTED  TRANSFORM  OF  FIRST  SWEEP 


phase  error  corresponds  quite  closely  to  the  broad  base  of  the  main 
lobe  of  the  transform.  For  comparison  purposes  the  weighted  transform 
of  an  ideal  rectangular  amplitude  and  linear  phase  signal  is  shown 
in  Figure  4.  The  time  delay  (range)  resolution  of  this  function  is 
4  [use c. 

Amplitude  and  phase  residuals  for  subsequent  sweeps  are  then 
computed.  With  the  first  sweep  acting  as  a  reference,  subsequent 
sweeps  are  compared  to  it  by  taking  the  ratio  of  the  amplitudes  and 
the  difference  in  the  phase  angles.  The  weighted  Fourier  transform 
is  then  computed: 


If  the  ionosphere  had  not  changed  its  character  during  the 

th 

interval  between  the  first  and  n  sweeps,  the  two  sets  of  data  would 
be  identical  and  the  transform  would  be  identical  to  that  of  Figure  4. 
For  the  data  being  discussed.  Equation  (4)  is  plotted  for  the  six 
subsequent  sweeps.  These  results,  shown  in  Figure  5  one  above  another, 
are  nearly  ideal. 

The  second  example  of  a  350  kHz  sweep  is  one  wherein  the  swept 
region  was  just  below  the  MUF  where  the  dispersion  is  quite  severe  and 
non-linear.  The  amplitude  fluctuation  in  Figure  6a  changes  frequency 
across  the  swept  region.  This  results  in  a  smearing  of  the  secondary 
lobe  after  transforming.  Note  also  how  the  amplitude  decreases  as  the 
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Figure  4.  TRANSFORM  OF  IDEAL  SIGNAL  WITH  40  db  TAYLOR  WEIGHTING 


Figure  5.  CORRECTED  TRANSFORMS  FOR  SWEEPS  2  THROUGH  7 
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Figure  6.  AMPLITUDE  ft  PHASE  RESIDUALS 


MUF  is  approached.  The  phase  residuals  to  a  first  order  fit  (Figure  6b) 
exhibit  the  same  variable  rate  fast  fluctuation  as  the  amplitude  data, 
as  would  be  expected  from  a  secondary  lobe.  In  addition,  the  more-or- 
less  triangular  shape  of  the  plot  indicates  that  the  delay  dispersion 
is  not  linear.  The  weighted  transform  of  the  second  sweep  is  shown 
in  Figure  7.  The  smeared  secondary  signal  and  broad  main  lobe  are 
particularly  evident.  The  transforms  of  Equation  (4),  which  simulate 
the  corrected  signals,  are  shown  in  Figure  8.  Even  after  10.5  seconds 
the  main  lobe  is  considerably  improved  while  the  secondary  signal 
has  not  been  affected  very  much  by  the  correction  process.  The  sweeps 
missing  were  discarded  due  to  equipment  malfunctions. 

Additional  data  were  recorded  at  Stanford,  California  on  June  20, 
1967  over  a  twenty-four  hour  period.  These  data  fall  into  two  cate¬ 
gories;  900  kHz  and  3  MHz  frequency  sweeps.  In  each  case  a  particular 
frequency  region  is  swept  repeatedly  for  20  seconds.  Of  a  total  of 
52  cases  of  900  kHz  data  two  examples  are  now  discussed. 

The  first  example  is  of  data  taken  about  9  FM  local  time  between 
10.5  and  11.4  MHz.  Amplitude  (A^)  and  phase  residuals  to  a  best 
straight  line  fit  (0^)  and  their  weighted  transform  are  shown  in 
Figure  9  for  the  first  sweep.  The  corrected  weighted  transforms  shown 
in  Figure  10  extend  for  18  seconds  after  the  origonal  sweep.  In  this 
case,  nearly  ideal  results  are  obtained  for  five  or  six  seconds  and 
quite  acceptable  results  for  the  full  18  seconds. 
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Figure  7.  WEIGHTED  TRANSFORM  OF  FIRST  SWEEP 
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Figure  8.  CORRECTED  TRANSFORMS  FOR  10  SECONDS  OF  DATA 
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Figure  10  CORRECTED  TRANSFORMS  OVER  18  SECONDS 


The  second  example  is  of  data  taken  about  9:30  AM  local  time 
between  16.0  and  16.9  MHz.  The  region  swept  was  right  at  the  MUF.  The 
amplitude  plot  in  the  upper  record  of  Figure  11  shows  severe  amplitude 
fluctuations  increasingly  close  together  as  frequency  increases  and 
then  tailing  off  as  the  MUF  is  passed.  The  dispersion  results  in  a 
1600  degree  peak  to  peak  phase  error.  The  transformed  signal  shown 
in  the  lower  record  of  Figure  11  is  spread  out  over  about  35  [asec. 
Corrected  transforms  one  and  two  seconds  later  are  considerably  better 
as  seen  in  Figure  12.  Further  corrections  10  to  12  seconds  after  the 
first  sweep  show  sidelobes  rising,  mostly  on  the  lefthand  side.  This 
unbalanced  result  is  typical  of  a  multi-path  situation.  The  missing 
sweeps  in  both  examples  are  due  to  equipment  difficulties. 

The  last  example  is  for  a  pair  of  3.0  MHz  sweeps  between  15.0 
and  18.0  MHz  taken  at  7:45  FM  local  time.  The  sweeps  were  5.0  seconds 
apart.  In  general,  the  effort  to  obtain  3  MHz  data  was  hindered  by 
the  limited  data  sampling  rate  imposed  by  the  maximum  stepping  rate 
of  the  digital  stepping  recorder  being  used.  The  receiver  was  limited 
to  a  bandwidth  of  50  Hz  which  was,  in  turn,  equivalent  to  a  50  |asec 
range  gate  due  to  the  time -frequency  mapping  of  the  linear  FM.  In 
this  case,  however,  the  receiver  output  stayed  within  the  50  Hz 
passband  over  the  3  MHz  sweep  region.  Amplitude  and  phase  between 
15  and  18  MHz  are  shown  in  Figure  13.  0^  ,  peak  to  peak,  is  5345  degrees. 
The  central  sections  of  the  uncorrected  and  corrected  transforms  are 
shown  in  Figure  14.  The  300  nanosecond  resolution,  inherent  in  the 
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lower  trace:  transform  of  upper  trace 
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Figure  12  CORRECTED  TRANSFORMS  UP  TO  12  SECONDS  LATER 


K> 

IO 

CVJ 

in 

<vi 

t 

<t 


o 

ro 

CJ 

in 

<\l 

i 

< 


UPPER  RECORD:  UNCORRECTED  FOURIER  TRANSFORM 

LOWER  RECORD:  CORRECTED  FOURIER  TRANSFORM  5  SECONDS  LATER 


waveform,  is  achieved  five  seconds  after  the  original  sweep.  These 
transforms  were  not  weighted,  thus  the  sidelobes  shown  on  the  correcteu 
transform  are  surprisingly  good. 

In  addition  to  the  processing  described  above,  which  basically 
consists  of  computing  Equation  (4)  for  each  case,  two  other  analysis 
techniques  were  employed.  Instead  of  comparing  the  first  sweep  with 
all  subsequent  sweeps  of  a  group,  each  sweep  was  compared  with  the 
next  one.  Thus,  the  following  weighted  transform  was  computed: 


This  analysis  technique  simulates  the  case  where  the  ionospheric 
transfer  function  is  measured  and  corrected  for  once  every  sweep.  It 
was  especially  useful  when  analyzing  a  disturbed  ionosphere.  A  resut. 
of  this  type  is  shown  in  Figure  15. 

In  another  test,  instead  of  a  best  straight  line  fit  to  the  pha^. 
data  a  third  order  fit  was  computed.  The  Fourier  transforms  using 
the  resulting  residuals  show  that  the  delay  dispersion  is  largely 
corrected  for  but  higher  order  fluctuations  caused  by  multipath 
phenomena  are  unaffected.  One  such  result  is  shown  in  Figure  16. 
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Figure  15  nTn  SWEEP  CORRECTED  BY  (  n  - 1  )th  SWEEP 
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Figure  16  TRANSFORM  OF  A,  AND  THIRD  ORDER  CURVE  rIT  ON 


SECTION  VI 


SUMMARY 

An  adaptive  signal -processing  technique  capable  of  compensating 
for  ionospheric  distortion  of  broadband  signals  has  been  described. 

A  computer  study  based  upon  data  gathered  in  a  Joint  MITRE-SEL 
experiment  with  the  SEL  coherent  sounder  has  verified  the  practical 
feasibility  of  the  processing  technique.  The  transfer  function  of  the 

ionosphere  was  measured  over  bandwidths  of  350  kHz,  900  kHz  and  3  MHz. 

\ 

These  data  displayed  both  dispersion  and  close  multipath  distortions. 
Application  of  the  corrective  processing  technique  to  these  data 
demonstrated  that  the  dispersion  can  be  corrected  for  periods  up  to 
eighteen  seconds.  The  more  rapidly  varying  distortions  associated  with 
the  close  multipath  (which  had  rates  up  to  one  cycle  per  second)  proved 
to  be  more  difficult  to  compensate  and  were  corrected  for  over  corres¬ 
pondingly  shorter  times. 
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